GRAVITATIONAL RADIATION FROM MONOPOLES 

CONNECTED BY STRINGS 

Xavier Martin and Alexander Vilenkin 
Institute of Cosmology, Department of Physics and Astronomy, Tufts University, Medford, MA 

02155, USA 

Monopole-antimonopole pairs connected by strings can be formed as 
topological defects in a sequence of cosmological phase transitions. Such 
hybrid defects typically decay early in the history of the universe but can 
still generate an observable background of gravitational waves. We study 
the spectrum of gravitational radiation from these objects both analytically 
and numerically, concentrating on the simplest case of an oscillating pair 
connected by a straight string. 



I. INTRODUCTION 

Monopoles connected by strings can be formed in a sequence of symmetry breaking 
phase transitions in the early universe The simplest sequence of this sort is 

G^HxU(l)^H. (1) 

For a semi-simple group G, the first of these phase transitions gives rise to monopoles 
which get connected by strings at the second phase transition. If both of these phase 
transitions occur during the radiation era, then the average monopole separation is al- 
ways smaller than the Hubble radius, and when monopole-antimonopole (MM) pairs get 
connected by strings and begin oscillating, they typically dissipate the bulk of their energy 
to friction in less than a Hubble time ||J3|]. 

A more interesting possibility arises in the context of inflationary scenario when 
monopoles are formed during inflation but are not completely inflated away. Strings can 
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either be formed later during inflation, or in the post-inflationary epoch. In this case, 
the strings connecting MM pairs can be very long. The correlation length of strings, £, 
can initially be much smaller than the average monopole separation, d; then the strings 
connecting monopoles have Brownian shapes. But in the course of the evolution, £ grows 
faster than d, due to small loop production, and to the damping force acting on the 
strings. Eventually, £ becomes comparable to the monopole separation, and we are left 
with MM pairs connected by more or less straight strings. At later times, the pairs oscillate 
and gradually lose their energy by gravitational radiation and by radiation of light gauge 
bosons (if the monopoles have unconfined magnetic charges). When the energy of a string 
connecting a pair is dissipated, the monopole and antimonopole annihilate into relativistic 
particles. 

The gravitational waves emitted by oscillating MM pairs add up to a stochastic 
background which can have an observable intensity @]. In order to calculate the spec- 
trum of this background, one first needs to find the radiation spectrum produced by an 
individual oscillating pair. This is our main goal in the present paper. 

In contrast to the case of cosmic string loops, the dynamics of MM pairs connected 
by strings has not been studied in any detail. We shall therefore concentrate on the 
simplest case of an oscillating pair connected by a straight string, for which the equations 
of motion can be solved exactly. Apart from its simplicity, this system has the advantage 
of being close to monopole-string configurations one would expect to find in the early 
universe. 

After reviewing the dynamics of monopoles connected by strings in the next Section, 



we calculate the gravitational radiation spectrum from an oscillating pair in Section [111 



then sum up and discuss our results in Section |TV|. Some technical details are given in 
the Appendices. 



2 



II. EQUATIONS OF MOTION 



The characteristic monopole radius 5 m and string thickness 5 S are determined pri- 
marily by the corresponding symmetry breaking scales, rj m and rj s . Typically, 5 m ~ t?" 1 
and 5 S ~ r/" 1 . In order to have monopoles connected by long strings, the two symmetry 
breaking scales should be well separated, i] s -C r) m , and thus the monopole radius is much 
smaller than the string thickness, 5 m <C 5 S . 

Assuming that the string length is much greater than its thickness, we can treat 
monopoles as point particles and strings as infinitely thin lines. The dynamics of a MM 
pair connected by a string can then be described by the action 

/ = — m J dsi — m J ds2 — /i / dS. (2) 

Here, \x is the string tension (which is equal to the mass of string per unit length), the 
first two integrations are over the monopole and antimonopole worldlines and the third 
is over the string worldsheet. 

The last term in Eq. (^) is the Goto-Nambu action for the string. Its variation gives 
the standard string equations of motion 

SIstring = ^ J <9 a ( V /II 77 at '^ b )5x ^Co^Cl ~ V J ^(a/ 3 ^ '^b Sx l*)dCo<Kl ( 3 ) 

where a = 0,1, ( a are a set of internal coordinates for the string worldsheet, x^ a = d a (x^), 
and '-fab = g iiv x il a x v h is the two dimensional worldsheet metric. The first term gives the 
equations of motion for the string which take a particularly simple form in the transverse 
traceless gauge 

i 2 + x' 2 = 1 , x • x = 0, (4) 

where the set of internal coordinates was taken to be it, a), and primes and dots refer to 
derivatives with respect to a and t, respectively. In this gauge, the dynamical equations 
are simply 

3 



x = x" (5) 

and can be solved exactly as 

x(*,<7) = i[a(<7-f)+b(<7 + *)]/2, (6) 

with the gauge conditions (|]) taking the form 

a' 2 = b' 2 = 1. (7) 

Since the string has two end points (one at each monopole), the spatial parameter must 
be restricted to an interval <7 2 (t)] where 

Xi(jt)=x(jt,(Ti(t)) (8) 

are the positions of the monopoles. 

The second integral in Eq. (]3|) can be turned into a boundary term and thus con- 
tributes to the variation of the monopole and antimonopole worldlines 

/ <9 a (V=77° 6 <6<^)<odCi = - / A la 7 a6 < 6 <^iM^i - / \2al ab x» h 5x 2 ,ds 2 , (9) 

where Ai a and \ 2a are unit vectors orthogonal to their respective worldlines and oriented 
into the string worldsheet, and x\ are evaluated on the monopoles worldline at (t, <Ji(t)). 
By definition, the same vector expressed in external coordinates is 

Kit) = A^V; = ± li {t)[a' i {t)x»{t,a i {t)) +*"*(*,*<(*))], (">) 

where 7j = (1 — x 2 ) -1 / 2 is the Lorentz factor of the monopoles. If we now add the terms 
coming from the variation of the monopole and antimonopole actions, we get the equations 
of motion for them in the form || 

A2 § - ^ <»> 
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where \x is the mass per unit length (and tension) of the string and m is the monopole 
mass. Since Af are unit vectors, it follows from flllD that 

a = fi/m, (12) 

is the proper acceleration of the monopoles. By multiplying equations (flT|) by x iu , it 
can be seen that only three of the four equations are independent. The time component 



equation of the system (11), expressing the exchange of energy between the length of 
string created or destroyed at the monopole end, and the kinetic energy of the monopole, 
takes the very simple integrable form 

a&i = ±7i, (13) 

while the spatial equations can be put in the form 

7 fx,(t) = ±a (j^j (t, a t {t)) = a^rv (14) 

Here, the Lorentz factor of the string at the location of the monopole, and 

rij is a unit vector pointing from the monopole in the direction of the string. 

The complete set of dynamical equations for the system of two monopoles connected 
by a string is thus given by the systems of equations (|5|) and fll4"|) with the constraint (^) . 
In general, the solutions of these equations are not periodic. 

Though the string part (Q-|5]) can be solved immediately in the standard form (0-0) • 
the motion of the monopoles is impossible to solve in general because of the presence of 
CTj(t), which is one of the unknowns, as a parameter of x' in equation (0). However, two 
sets of particular exact solutions can be found by assuming either that a(t) is constant or 
that x' does not depend on its spatial parameter, thereby removing the problem. 

In the first solution, the string has the form of a rotating rod of length / with the 
centrifugal force acting on the monopoles balanced by the string tension: 
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x(t, a) = Rsin(a/R)e(t/R), (15a) 
Xi(t)=±(l/2)e(t/R), (15b) 
<7i(t) = ±i?arcsin(//2 J R), (15c) 

where e(9) = (cos 9, sin 9) is the radial unit vector associated with the angle 9 and 
/ = ((1 + 4a 2 -?? 2 ) 1 / 2 — l)/a. We note that in the case when the monopoles are rela- 
tivistic, that is when aR 3> 1, the energy of the string is much larger than the energy 
of the monopoles, E s /E m ~ aR. We will make use of this property to approximate the 
gravitational radiation spectrum of this solution by that of a simple cosmic string loop 
(rotating double line). 

The second solution describes an oscillating pair of monopoles connected by a 
straight string: 

x(£, a) = ae, (16a) 
x .( t ) = ±^^(7o - v/l + (7o^o-^|) 2 )e (16b) 
ai(t) = Xi(t). (16c) 

Here e is a unit vector along the string, which we choose to be directed along the x-axis, a is 
the monopole proper acceleration defined in (|T^), t> and 70 = (1 — Vq)~ 1 / 2 are respectively 
the maximum velocity and Lorentz factor of the monopoles, reached at t = 0. The 
monopole and antimonopole meet at t — and could be expected to annihilate. However, 
this solution is considered as an approximation for an almost straight configuration where 
the monopole and antimonopole would merely come close to each other and would not 
collide. Besides, as we already mentioned, the monopole radii are much smaller than the 
string thickness, thus the monopoles are not likely to collide, even for a straight string. A 
peculiar feature of the solution (|l^) is that the monopole accelerations abruptly change 
direction when the monopoles meet and pass one another. 

The solution (^) is valid for \t\ < 70^0/0. At t = —•yoVo/a, the monopoles are at 
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rest, with the string having its maximum length, L = 2(70 — I) /a. At t = +7ot>o/a, the 
monopoles come to rest again, with their positions interchanged. As far as gravitational 
effects are concerned, since the monopole and antimonopole have the same mass, they 
are identical, and the full period of motion is T = 2j vo/a. On the other hand, if 
electromagnetic effects are considered, the monopole and antimonopole have opposite 
charges, and thus equation (|T6|) describes only half a period, T = 47 t> /a, the other half 
period being obtained by exchanging the positions of the monopole and antimonopole: 

x 1 (t + T/2)=x 2 (t) = -x 1 (t). (17) 

In the following section, we shall study the gravitational radiation from an oscillating 
pair described by the solution ([16]). Radiation from the rotating rod configuration will be 
discussed in Appendix |A]. 

III. GRAVITATIONAL RADIATION 

The power in gravitational radiation from a weak, isolated, periodic source to lowest 
order in G, can be found from the following equations, without any further assumptions 
about the source 

P = Y,Pn = Y,[d^, (18a) 

HP G } 2 1 

= ^(T;> n ,k)T^V„,k) - -|T/K,k)| 2 ). (18b) 

Here, dP n /dQ is the radiation power at frequency u n = lixnjT per unit solid angle in the 
direction of k, |k| = u n , T is the period of the oscillation, and 

T^(cj„,k) = i [ T dtexp(iu n t) [ d 3 x exp(-ik ■ x)T^(x, t) (19) 
T Jo J 

is the Fourier transform of the energy-momentum tensor. 



For the solution (|i"6[ ) considered in this section, since the system is one dimensional, 
the energy-momentum tensor has only three non zero components: T 00 , T 01 and T 11 . It 
also satisfies the conservation equations V M T MW = which in Fourier space can be written 
simply 

u n T 0u = k l T iv . (20) 

This means that T^ v has in fact only one independent component; the simplest choice is 

T 01 (t,x) =m^ v -a\t\){5(x-x 1 (t)e)-5(x + x 1 (t)e)], (21) 

which has no contribution from the string part of the system. Its Fourier transform as 
defined by fli~9|) can be simplified to 



T 01 (uj n , k) = m^v I n {u), (22a) 
I n {u) = f 1 £d£[cos(n7r(l - e - - + u^e + V(7o^o) 2 )) (22b) 

J v v 

- cos(n7r(l - e + ^ - u^e + l/(7o^o) 2 ))], 

where cu n = ri7ra/ (70^0) is the angular frequency of the n-th mode and we have introduced 
the notation u = k x /u n . The two other non-zero components of the Fourier transform of 
the energy-momentum tensor can then be deduced from the conservation equations (|20| ) 

as 

T°V n ,k) =uT 01 (u n ,k), (23a) 

T n K,k) = -T V„,k). (23b) 
u 



The gravitational energy radiated in the mode n can then be expressed from ( |18b|) as 



P n = 2G(mrfi) 2 C du(- - uf\I n {u)\\ (24) 
Jo u 



where I n (u) is given by (|22aj) and depends only on u = k x /uj n . This power spectrum can 



not be integrated in a closed form. However, it is possible to get analytic approximations 
at low and high frequency. 



At high frequency, an expansion of (22) in 1/n can be made as shown in Appendix 
First, a change of variable 



C = ±£-W£ 2 + 1 /(7cA)) 2 (25) 



enables us to rewrite (E2b|) in the standard form (iBll 



W = 7j ^ / ° C0S Ml + - + OR, (26a) 



(1-U 2 ) 2 J-i-jL ' V 

1 — u 2 -, /0 1 — u 2 , ,_<, 1 — u 2 



f(u, C) = 2(C 2 + -^rf 12 ~ VHC 2 + -^r)- 112 - (26b) 

7o^o 7o^o 7o^o 

The expansion is obtained by repeatedly integrating by parts in Eq. ( |26a|) . The dominant 
term is 

/ ( U ) „ uv ° ( _J L^l l-v 2 o ) m 

n{) (1 - u 2 ) \1 - u 2 v 2 (l + uv Q ) 3 (l-^o)V' 

To get the interval of validity of this expansion, this term must be compared with the 
next as shown in Eq. ( [Bo] ). For simplicity and because it is the most interesting case, 
we shall assume in the following that the monopoles are ultra-relativistic, 70 ^> 1. Then, 
gravitational radiation is beamed in the monopole's direction of motion, into a cone of a 
small opening angle 9 ~ l/7o- Thus, the main contribution to the total power fll8|) comes 
from values around u ~ 1 and the comparison of the two terms can be performed there. 
This gives 

(1 - « 2 )7$, (28a) 
If'CS ^(l-« 2 )7o 8 , (28b) 

VQ 

so that the expansion is valid as long as 

n » 7o- (29) 
The dominant term of the power spectrum can be expressed as 



p 2Gu 2 /-i 
n z n z vfi Jx-vo 



II : I I I / / 

(30) 



Avl %w - 1 1 + v 2 - w 

< AO' 



w(2 — w) 7q^ 3 7q(2 — w) 3 

The 1/n 2 behavior at infinity could be expected since the second derivative of the stress 
energy tensor has a discontinuity at t — 0. When 7 is also assumed to be large, that is 
the monopoles reach ultrarelativistic velocities, the integral (|30|) can be simplified to 

P„«|i(^)V (31) 

For n 7q with j 3> 1, it is legitimate to neglect the terms (70^0) ~ 2 m the square 
roots of equation ( |22b| ). This gives 

r^m , 1 \ /sin(n7r(l — u)/2) \ 2 

T^ n ,k «8m 7 o« — Si £-L\ , 32 

y 7277(1 — U Z J J 

and the power in the n-th mode can be simplified to 

P n « — G/i / du— — . 33 

nir Jo u z (l — — r 

It is interesting to note that this low frequency behavior is independent of the maximum 
Lorentz factor of the monopole 70, as long as it remains large. Furthermore, for large n, 
it is possible to neglect ii/titc and extend the integration to infinity so that we get 

16 2 /*°° sin 4 u 4Gu 2 
P n « — Gu 2 / du — — = — — . 34 

nix Jo u z n 

The low frequency behavior (|3^) has been plotted in Figure [1] and is indeed well approxi- 
mated by 4Gu 2 /n for n ^ 30. From a cosmological point of view, the quantity of interest 
is the gravitational energy per logarithmic interval nP n which is therefore quasi-constant 
in the frequency interval 7q > n ^ 30. Curiously, the low frequency spectrum for the 
other solution, the rotating rod fll5D , though it is very different from the straight string 
solution, also behaves like 1/n (see Appendix [A]). 

The power spectrum can be computed numerically by integrating successively ([22] ) 
and (HU). The evaluation was done for various values of 70 and mostly exhibits a smooth 
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evolution from the low frequency (|33|) to the high frequency (^T|) behaviors. An example 
with 7q = 25 is shown in Figure [l]. 

It is also useful to know the behavior of the total gravitational energy loss rate P 
as a function of 70. Once again only a numerical solution is possible. However, instead of 
adding up the power in different modes, it is much faster to compute P directly without 
going to Fourier space. The calculation is outlined in Appendix [C], and the resulting 
radiation rate is plotted in Figure |2| as a function of 70 . Empirically, it can be closely 
approximated by 

P( 7o ) ^ (81n( 7 o) + 2.2)G> 2 , (35) 

which is consistent with a spectrum behaving first like 4Gn 2 /n up to n ~ 7q and then like 
(64/57r 2 )(7o/n) 2 G/i 2 . For large values of 70, the flat low frequency part of the spectrum 
makes the dominant contribution to the total gravitational power emitted. In the case 
70 = 25 considered above, the power computed numerically is P = 28.0G/1 2 (in agreement 
with a summation of all the modes of the spectrum found in Figure P while the algebraic 
approximation gives P ~ 28.3G/i 2 with a relative error of only 1%. 

IV. CONCLUSIONS 

We have analyzed the gravitational radiation from an oscillating monopole- 
antimonopole pair connected by a straight string. The motion of the pair is described by 
Eq. flTED. The gravitational radiation is emitted at a discrete set of frequencies uj n = nuj\, 
where lo\ = Tta/joVo, a = fi/m is the proper acceleration of the monopoles, /x is the string 
tension, m is the monopole mass, vq is the highest velocity reached by the monopoles, and 
70 = (1 — f 2 ) -1 / 2 is the corresponding Lorentz factor. In the most interesting case of ultra- 
relativistic motion, 70 3> 1, most of the radiation is emitted in the range 1 ^ n 7o> with 
a spectrum P n m 4Gfj, 2 /n. For n>l, the spectrum can be approximated as continuous 
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with 

dP/duj « AGfx 2 /cu. (36) 

At higher frequencies, dP/du oc u~ 2 . The total radiation power is 

P^8G/i 2 ln( 7o ). (37) 

The one dimensional solution ( |I6|) should be regarded as an approximation for a 
more general configuration of monopoles connected by a nearly straight string. In a more 
realistic case, the 1/uj and 1/uj 2 behavior is expected to be modified for sufficiently large 
n ^> n c .The characteristic value n c and the corresponding frequency uj c can be estimated 
from uj c ~ (A/) -1 , n c ~ //A/, where A/ is the monopole separation at which deviations 
from the straight-line shape become important and / is the maximum extent of the string. 
Typically, Al is comparable to the minimum distance between the monopoles as they pass 
one another. 

We cannot tell from our analysis how the spectrum is modified at uj > uj c . This 
remains a problem for future research. We expect that eventually P n will fall off exponen- 
tially at n — > oo, but there may also be some intermediate regime. Since no solutions of 
the equations of motion are known in which the string would deviate from a straight-line 
shape, this problem will probably have to be tackled numerically. In particular, one could 
employ the numerical simulations of the monopole-string system that are currently being 
developped ||. 

Though in the simplest models of symmetry breaking all the magnetic flux of the 
monopoles is confined in the strings, in more realistic models, stable monopoles can have 
unconfined magnetic charges. In such a case, monopoles can lose energy by radiating gauge 
quanta. The gauge fields associated with the magnetic charge may include electromagnetic 
or gluon fields, but may also correspond to broken gauge symmetries and have non zero 
masses. The gauge boson radiation is important since it can greatly affect the lifetime of 
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the pair and thus the total output of gravitational waves. An evaluation of the radiation 
of massless or massive gauge bosons by monopoles will be given elsewhere ||. 
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APPENDIX A: GRAVITATIONAL RADIATION SPECTRUM FOR A 
"ROTATING ROD" SOLUTION 



The rotating rod solution ( Tj5|) does not seem likely to arise naturally in the early 
universe but could give some insight into what happens for a configuration very differ- 
ent from the one- dimensional solution ([161) . We shall concentrate on the case when the 
monopoles are relativistic, that is when aR ^> 1. In that case, it is easy to check that the 
energy of the string is much larger than the energy of the monopoles, E s /E m ~ aR ^> 1. 
This means that the monopoles can be ignored in calculating the gravitational radiation 
of the system. As for the string part of the system, it should be well approximated by 
the straight rotating double line solution 

x(t,cr) = Rsm(^)e(t/R), (Al) 
R 

with \a\ < ttR, as long as the wavelengths of the gravitational waves remain large com- 
pared to the length difference of the rotating string segments in the solution (|15|) and its 



approximation n <C aR. Since the loop ( |A1[ ) is made of two straight strings, its 

mass per unit length \x should be half that of the straight rod solution: 

pL = n/2. (A2) 

The straight loop solution ( |A~l|) is singular at its "end points" a = ±ttR/2, which move 
at the speed of light. This singularity results in a 1/n decay of its gravitational radiation 
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spectrum and thus in a divergent total power. This is not really a problem since our 
approximation is only valid at low frequencies, n -C aR. However, it is for this reason 
that, though the gravitational radiation of similar loop solutions has been studied in the 



past the calculation for this particular solution do not appear to have ever been 

done. 

Using the standard formula for gravitational radiation from a string in the direction 



of the unit vector k, [12 



dP 

^(k) = 87rG/iV{|/ n (n 1 ) J„(m) - J n (n 2 ) J n (n 2 )| 2 + (A3) 



|4(ni)J n (n 2 ) + 4(n 2 )J„(ni)| 2 } 



where rix and n 2 are unit vectors orthogonal to k and to one another, J n (n) and J n (n) 
are defined as 

1 r*- 1 nidi 
J B (n) = - / dCa'CC) " nexp[- — (C + k • a(C))], (A4a) 
L Jo 2 

J„(n) = l / L rfCb'(C)-nexp[^(C-k-b(C))], (A4b) 
L Jo 2 

a(£) and b(£) are the generators of the string worldsheet as defined in equation (|6|), L is 
their period and u is the angular frequency of the loop. For the loop solution (KI), after 



integrating the gravitational power over all directions and replacing \x by its value (A~2), 
we have 

n q o / ^ COS OC -j. a t COS OC / k \ 

P n = 87r G/i n I smxdxl — j — J (nsini) +6 — ~ — (A5) 
jo sin x sin x 

J 2 (n sin x) J n ' 2 (n sin x) + sin x)J n ' 2 {n sin x)} } 

where J n {x) are Bessel functions of the first kind. This integral can be computed nu- 
merically and gives a power spectrum which behaves almost exactly like 5.75G/i 2 /n, even 
at low frequencies. The spectrum of the rotating rod solution ( |I~5l) should therefore be- 
have like 1/n for low-frequency modes n aR, but is expected to fall exponentially at 
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high frequencies (since the stress energy tensor of the system, and all its derivatives, are 
regular). 



APPENDIX B: LARGE-n EXPANSION 

In this Appendix, we develop a method for deriving large-n asymptotic expansions 
for integrals of the form 

/•n.4-2 

f(u) cos(n7r(-u — a))du, (Bl) 

where / is a regular function which does not depend on n. Examples of such expressions 
are found in this paper in equation (^6|) , and more generally in problems involving Fourier 
transforms. The trick is to integrate by parts the cosinus 

1 1 r a + 2 
I n = — [f{u) sm(rnr(u — a))]" + / f'(u) sin(n7r(w — a))du 

T17T 717T J a 

1 ra+2 

[f{u)cos{rm{u-a))] a a +2 - f"(u) cos(nn(u - a))du. (B2) 



n 2 ir 2 



The last integral is of the same general form as the initial expression (|B1| ) so that this 
expansion can be iterated to get 



In = Uf^ia + 2) - / (2Z - 1} (a))^r- (B3) 

In this paper, we are only concerned with an expansion for large n to the first significant 
order 



[f^Wr + Oi—-). (B4) 



The interval of validity of this expansion can be evaluated by requiring that the first 
neglected term in the expansion be much smaller than the one we keep. This gives 

" 2>> W- (B5) 
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APPENDIX C: DIRECT CALCULATION OF THE TOTAL GRAVITATIONAL 

POWER 



Here, we derive an expression for the total gravitational power radiated by the one- 
dimensional solution flUf). The metric perturbation W = g^ v — rf" induced by a weak 



source at a large distance r in the direction of the unit vector n can be expressed as 

AG 



Air f 

h» u (t, rn) = — / T»"(t r , x')dV, (Cla 

hT = - ^h p p rf v , (Clb 

t = t r + |x' -x(t r )|, (Clc 



where ff is the Minkowski metric. The gravitational power emitted by the source can 
then be expressed as 



dP , s r 2 /T .„- 1,- 

— (t) = 

dQ K J 32nG 



(^^V-rfW. (C2) 



which no longer depends on the distance to the source r. In the present case, there are 
only three non zero terms h pu given by 

h ™ = ho + (^k + (C3a) 

2 y du du J 

h* = ™(^ + *£), (C3b) 
2 y du du J 

.11 h mn x fdul du\ l + 2u 2 2r 1 + 2u 2 lr \ 

h = hi — — 1 — ) + ml 1 , (C3cJ 

2 y du du 72 — n x U2 r li + n x u\ r J 

where ho and hi are constants which do not contribute to the total gravitational power, 
uir — lo v o ~ d\ti r \ and U2r = 7o^o — a\t2r\ are functions of the retarded time for each 
monopole which can be expressed as functions of u = 70 Vq — a\t\ by 



x\ - Vi 



(-iyr^yJvf + 1-nl 
1-nl 

(-1)^70 - u. (C3e) 



u ir = y — - , (C3d) 
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What we are really interested in is the gravitational radiation power averaged over a pe- 
riod. This is obtained by integrating equation (|C2|) over all directions and then averaging 
over time. This gives 

r 2 novo ri ^oo ^11 -0i ^oo _l01 

P = — / du dnJh +h +h )(h +h -h ). (C3f) 

16G7 fo Jo Jo 

Though the expressions ( |U3| ) are complicated, finding the average total emitted power 
P(7o) only requires a double numerical integration. Since the evaluation of each point of 
the spectrum also required a double integration, this method is much more efficient than 
simply summing the power of all the modes. 
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FIG. 1. A log-log plot of the gravitational radiation spectrum nP n /(Gfi 2 ) in the case 70 = 25 
(solid line) with its high frequency approximation (31) (dashed line) and low frequency approx- 
imation (|33|) (dotted line). 



FIG. 2. Total gravitational power emitted (solid line) and its empirical algebraic approxima- 
tion (35) (dashed line) as functions of In (70). 



18 



Figure 1 



4.0 



3.0 



2.0 



_ 1.0 

CM 

o 

Q_ 

^ 0.0 



-1.0 



-2.0 



-3.0 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 

ln(n) 



